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Abstract 

We discuss the collective dynamics of self-propelled particles with selective attraction and repul- 
sion interactions. Each particle, or individual, may respond differently to its neighbors depending 
on the sign of their relative velocity. Thus, it is able to distinguish approaching (coming closer) and 
moving away individuals. This differentiation of the social response is motivated by the response to 
looming visual stimuli and may be seen as a generalization of the previously proposed, biologically 
motivated, escape and pursuit interactions. The model can account for different types of behavior 
such as pure attraction, pure repulsion, or escape and pursuit depending on the values (signs) of the 
different response strengths, and provides, in the light of recent experimental results, an interesting 
alternative to previously proposed models of collective motion with an explicit velocity-alignment 
interaction. We show the onset of large scale collective motion in a subregion of the parameter 



x. 

space, which corresponds to an effective escape and/or pursuit response. Furthermore, we discuss 
the observed spatial patterns and show how kinetic description of the dynamics can be derived 
from the individual based model. 
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I. INTRODUCTION 



Collective motion in biology, as observed for example in flocks of birds, schools of fish or 
within bacterial colonies, is a fascinating display of natural self-organization. Over the years, 
it has been the topic of numerous scientific publications addressing it from very different 
angles and with different questions in mind, both experimentally and theoretically. From 
a more biological perspective the interesting questions are the evolutionary advantages and 
the biological and ecological function of collective behavior in various species jl-4], whereas 
physicist addresses rather the question about universal laws and phase-transition behavior 
in minimal models of collective motion js- 12] . The design, control and stability of collective 



dynamics in multi-agent systems is also a major research topic in engineering 



13 



-]l6|, and 



the general properties of related mathematical models are under active investigation in 



mathematics 



13, 



18| 



Most of the mathematical models for collective motion proposed in the literature, contain 
some sort of explicit velocity-alignment mechanisms, which tends to align the velocity of a 
focal individual with the velocity of its neighbors [l, 5, 



19 



22|. However, recent experimental 



studies of collective behavior in fish do not find any clear evidence for the existence of an 



explicit velocity-alignment interaction 23|, ]24j . Only relatively few models, have analyzed 
the onset of collective motion without an explicit alig nment mechanisms based on purely 



repulsive and attractive interactions (see e.g. 25N29|). Recently, motivated by empirical 
evidence for cannibalism as the driving force of collective migration in certain insect species 



30 
4 



3l| , we have proposed a model of collective motion based on escape and pursuit responses 



11| . In this escape-pursuit model, individuals are reacting to their neighbors by moving 



away from others approaching them from behind (escape), and/or increasing their velocity 
towards those who are moving away in front of them (pursuit). This kind of social response 
requires individuals to distinguish between approach and movement away as well as between 
individuals in front and behind them. This previous model can be considered to belong to a 
broader class of selective attraction-repulsion models, which we believe are very promising 
for theoretical modelling of collective motion in biology. Here we discuss and analyze a 
generalization of the original escape-pursuit model to the case where self-propelled agents 
(or particles) are responding selectively to approaching and moving individuals without 



taking their relative position into account 



32] . Furthermore, we do not put any restrictions 



on the sign (direction) of the effective social forces modelling the selective response. This 
allows to account for different social behavior types, such as pure attraction, pure repulsion 
and escape and pursuit, in a single model, by the same set of social forces, only by changing 
the values of the response strengths. 

We will start with the definition of the individual based model in terms of stochastic 
differential equations and continue with the derivation of a kinetic description of the system. 
Finally, we will discuss simulation results with a particular focus on the emergence of large- 
scale collective motion. 

II. INDIVIDUAL BASED MODEL 

We consider a system of N self-propelled particles in two spatial dimensions, which move 
with a constant speed sq in a spatial domain of size Lx L with periodic boundary conditions. 
The interaction between different particles (individuals) is modelled as an effective social 
force Fj. The evolution of the system is determined by the following equations of motion 
for the positions rj and the polar orientation angle ip, which determines the direction of the 
heading unit vector e h;i (t): 



The temporal evolution of tfi is determined by the turning of the individual due to social 
interactions F iytp and random (angular) fluctuations with the intensity D v . The (angular) 
social force is given by the projection of the total social force vector Fi tV> = F^e^ on the 
angular degree of freedom with e Vi = (— sinyjj, cosy2j). The angular noise £ v is Gaussian 
white noise with zero mean and vanishing temporal correlations. 

The total social force is given by a sum of three components: 



The first term represents a short range repulsion responsible for collision avoidance. It reads 




(1) 



(2) 



F 




(3) 




(4) 
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with fi r > being a constant repulsive turning rate. The Heaviside function 9{l T — r^) 
ensures that the repulsion takes place only if the distance = |rj — r»| between the focal 
individual i and the respective neighbor j is below the repulsion distance l r . The total 
repulsive response is normalized by the number of individuals within the repulsion distance 

JV 

Nr = Nr{t) = Y^6{l r -r 3l ). (5) 

i=l 

The other two forces read: 

1 N 

f m = Jf-ny^PalVjilTjidils ~ T^Ofai ~ l r )0(+Vji), (6) 
1 N 

f a = ^"^y ^2 Vrn\Vji\r jiO{l s ~ T^fai - l r )6{-Vji). (7) 

Both forces represent averaged two- individual interactions, which act always along the unit 
vector pointing towards the center of mass of the neighboring particle = (ij — Ti)/\Tj — rj|. 
The first one, f m , represents the response to approaching individuals characterized by a 
negative relative velocity Vji = (v, — Vj)fjj < 0. The second, f m , is the corresponding 
response to moving away (or receding) individuals characterized by positive relative velocity 
Vji > 0. This differentiation is reflected by the last Heaviside functions 6(±Vji). The two 
other step functions are identical for both interactions and restrict these social responses 
to neighbors within a sensory range l s but outside the repulsion zone. The parameters 
/i m>a determine the turning rates due to the respective interaction. Both force terms are 
proportional to the relative velocity, which lead to stronger responses to faster approaching 
or receding individuals. Furthermore, they are normalized by the respective number of 
individuals for the corresponding interaction type: 

N N 

N m {t) = Y,0{r jl -lr)e{l s -r ]l )9{+v ]l ), N a (t) = ^ 6{ rji - l r )6{l, - r^-v^. (8) 

i=l i=l 

Here, we used for simplicity step-like functions for the spatial dependence of the different 
interaction. The general results will not be altered by other smooth functions of the distance 
as long as they decay sufficiently fast in order to ensure local interactions. Please note that 
the definition of the step-like interaction zones resembles the two-zone model introduced by 
Couzin and co-workers 



331 ] . However, the model discussed here does not contain an 



explicit velocity- alignment. 
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FIG. 1. Schematic visualization of social interactions: The focal individual i (gray) can interact 
with individuals within its sensory range l s . Hereby it distinguishes between "approaching" (red) 
and "moving away" (blue) individuals. The decisive factor in the distinction is the sign of relative 
velocity Vji defined by the projection of the velocity difference of neighbor j and the focal individual 
(vjj = Vj — Vj) on the relative position unit vector rji = Tji/\rji\. 

A schematic visualization of the interaction scheme with the differentiation between ap- 
proach and moving away is given in Fig [TJ 

The social forces f a / m can lead independently to a repulsive (attractive) response to 
approaching individuals for \i a < (/i a > 0) and a repulsion (attraction) to individuals 
moving away /i m < (/i m > 0). In the /i m/ u a -parameter space we distinguish the four 
quadrants corresponding to different behavior types (see also Fig. [2]): 

i. Pure Repulsion: repulsion from approaching and moving away individuals: fi a < 
and /i m < 0. 

ii. Escape and Pursuit: repulsion from approaching individuals fi a < 0, attraction to 
moving away individuals /i m > 0. 

iii. "Head on Head": attraction to approaching individuals fi a > 0, repulsion from 
moving away individuals /i m < 0. 

iv. Pure Attraction: attraction to approaching and moving away individuals: fi a > 
and /i m > 0. 
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There exist also the special cases with fi m / a = and fi a /m > (/x a /m < 0), which correspond 
to a selective attraction (repulsion) only to approaching/moving away individuals, and the 
case of particles interacting only via short-range repulsion (/x m = fi a = 0). 

We refer to the situation fi a < and fi m > as "Escape and Pursuit", due similar 



behavior as in the original Brownian particle model [ll|. For /x a > and \i m < the social 
forces lead to a preference to move towards other individuals which already are coming closer 
and therefore favor (in particular at low densities) frontal collisions between individuals. We 
refer to this regime as "Head on Head" . 

III. KINETIC DESCRIPTION 

In this section we derive a kinetic description for the above individual based model. For 
this purpose we introduce the iV-particle probability density function (PDF) 

Pn(ti, <fi, r 2 , tfi2',---\ ?n, ¥n] t) , 

which determines the probability to find a particle (individual) at time t, at position r i; with 
velocity pointing in direction tpi (i — 1, 2, . . . , N). It is normalized with respect to integration 
over space and over all angles. Further on, for simplicity, we assume that correlations 
between particles can be neglected. Therefore, the iV-particle distribution density shall 
factorize, i.e. Pn = Ilf =l P(ri,(pi,t). In agreement with flTJ and (j2J), we can write down the 
Fokker-Planck equation (FPE) for the PDF of the z-th particle 

with Bfr. = (cos ifii, simp i) T being the unit vector in the heading direction of individual i, 
and e Vi = (— sin<^, cos ipi) T being the angular unit vector perpendicular to e^. The above 
FPE is nonlinear since the interaction force Fj depends on the probability density for the 
position and the velocity angle of the particles within their sensory range. 

We will now reduce the description to moments of the one-particle PDF, which are the 
particle density p(r, t) 

/*27T p 

p(r,t) = drP(r,<p,t), /drp(r,t) = l, (10) 
Jo Jv 

and the expectation values of the cosine and sine of the velocity angle defined as 

/*2"7r /*27r 

s(r u t) = / d<psm((p)P(ip,t\ri) , c(r u t) = / dip cos(<p)P(y?, t|rj) . (11) 

Jo Jo 
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FIG. 2. Examples of spatial configurations for different regimes: (A) "Escape and Pursuit" (fi a = 
-3.0, fi m = 3.0), (B) "Pure Attraction" (/x = 3.0, fi m = 3.0), (C) "Pure Repulsion" (//„ = -3.0, 
= —3.0) and (D) "Head on Head" (/i a = +3.0, /i m = —3.0). The arrows and their color indicate 
the direction of motion of individual particles. The mapping of the color to the directions is shown 
as an inset in (A). The different panels are arranged according to the location of the corresponding 
regime in the interaction parameter space, with the origin (/x m = [i a = 0) being in the center. 
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where the conditional probability density function of the velocity direction ip is defined 
through the relation P(r,tp,t) = P((p,t\ri)p(r,t). Similar approach was previously used in 



the context of swarming of Active Brownian Particles in [34j-|36]. 

Integrating the Fokker-Planck equation over the velocity-angle yields the continuity 

equation which reads 

d d d 

— p{vi,t) = - s Q —c{v h t) p(r h t) - sq— s(r h t) p(r h t) . (12) 

Similarly, we derive the equations for the angular moments 

d d d 

— c{r h t)p{r h t) = - s — ■c 2 {r i ,t) p{r h t) - s —cs(r h t) p{r h t) 

-—%-c(r h t)p{vi,t) - — (Fie^simpi) p(ri,t) (13) 
So s 

d d d 

— s(r h t)p(r i ,t) = -Sog^-cs{ri,t)p(ri,t) - s ^— s 2 (ri, t) p(v h t) 

- : %c(r i ,t)p(r i ,t) + — (Fie^ cos ipi) pfat) (14) 
s o s 

Here, we introduced the averaging over the conditional PDF as (...) = J Q 27T d<p ... P(<p, t\r) 
which are functions of the position r and time. The second moments of the trigonometric 
functions of the direction angle read 

c 2 (r;,t) = (cos 2 (<^)) = -(1 + (cos(2pi))) « -, 
8 2 {T i} t)= (sin 2 (^)) = -(l - (cos(2^)» « ^> 

cs(Ti,t) = (cos((fi) sin((fi)) = - (sin(2<pi)) « 0. (15) 

The approximation for the second moments indicated on the r.h.s. of the above equation, 

will be used later on. It neglects higher harmonics in the equation of motion for the angular 

PDF, which essentially restricts the analysis to the relaxation of the slowest (fundamental) 

modes of the angular PDF. 

In order to obtain an equation for the density, we insert the equations of the angular 

moments into the continuity equation. For this purpose, we take the temporal derivative of 

the spatial density for a second time 
d 2 D d \ 



d 2 d 2 d 2 

■c 2 (r i ,t)p(r i ,t) + 2 cs(Ti,t)p(Ti,t) + — I s 2 (r i ,t)p(r i ,t) 



oxf oxidyi oyf 

1 d . . . . 1 d . . . . 

H 7^— (Fje^ sm y^) pfa, t) — (F^ cos y3i)p(r;, t) 

so oxi so dyi 



If we introduce in this equations abbreviations for the variances of the second angular mo- 
ments defined as 



T XhXi (ri,t) = 02(1-3, t) - c 2 (r u t) , 
T XiiVi (ri,t) = csfat) - c(r i5 t) s(r i5 t) , 

^w.w(ri,<) = s 2 (ri,t) - s 2 (r;,t) , (17) 



the equation for the density becomes 
d 2 D,„ d 



<£_ 

dt 2 ' S 2 dt 



Id Id 

H ^— (Fie w sin p(r i? t) — (F^ cos ^Pi)p{r u t) . (18) 

So OXi Sq oyi 

We note, one still needs equations or expressions for the variances T u ._ v . which have to be 

inserted. 

Similarly, we separate the density from the angular moments using the continuity equa- 
tion. After a few steps, one obtains the equation of motion for the mean cosine 



d . . d , . d \ , . (Fie w sin<£i) D v 



SO 



and for the mean sine 



f d d 

i^T^(Ti,t) p{Ti,t) + —T Xi} y.{v h t)p{v h t) 



(19) 



•9 , ^ d i ^ 9 I ^ i (Fie w cos (p^ D v 



dt dxi dyi J s sg 



(20) 

respectively. We multiply the first equation f lT9|) by the mean cosine and the second one 
( 120]) by the mean sine and obtain the equation for the order parameter: 

i|^ + C (r,,t)A + s ( r4; t)^_| ( c 2 (r 4 ,t) + s 2 (r 4 ,t)) = - ^(c 2 (r 4 , t) + s 2 (r 4 , t)) 
- — (c(ri,t) (sin(v9 i )F i e ¥ , ! ) - s(r i5 1) (cos(v? i )F i e (Pi ) ) 

" 7^{£ T -<«>«>><«>«> +|"^, w ( r ,i)p( r „i)} (21) 
sps(rj,t) \ d 

The only term, which can induce an instability of the homogeneous, disordered solution 
in Eqs. f lT6|) and (l2Tj) contains the interaction force. The remaining terms, which contain 



f d d 1 



9 



only the zeroth up to second moments of the orientation, are present also in the case of 
non-interacting particles. They describe the relaxation dynamics of the resulting patterns 
towards a steady state, but are not the source of an eventual inhomogeneous solutions. 

We see that a further analysis needs still the treatment of integrals which contain the 
interaction forces, which can be very tedious in the general case. Further below, we will 
consider principal cases where \x m and \i a have the same absolute value and differ only in 
their sign and derive approximations for the respective integrals. 

The force changing the direction of the ith-particle can be formulated by introducing 
an interaction parameter depending of the distance between particles and their relative 
velocities. We define 

{-§• for < < l r 
(22) 
fj, z (vji) for l r < rji < l s 

where is the distance between from the zth to the jth particle and Vji again the relative 

velocity projected on the distance vector. Otherwise, this function vanishes. The « z have 

to be specified for positive and negative relative velocities. They read, in agreement with 

© and ©, 

= { Na ■ (23) 

|+M for%>0 

The normalizing numbers in denominators can be expressed as integrals over the respective 
distances and relative velocities. In detail, with the definition of the total number density 
as uq = N/V, where V = L 2 is the total volume (area in 2d) one gets 

p /*27r 

/ dYji j dtpjP(<pj,t\r + Tji) p(r + r ji , t) , (24) 
Jav z Jo 

where the AV z ,z = r,a,m are the (effective) partial volumes for the different interaction 
types, respectively. For the integration one can transform to polar coordinates, integrating 
over the distance and the polar angle 0£j, determining the position of particle j with 
respect to the focal particle i. 

The force in the third term on the r.h.s. of the Fokker-Planck equation (FPE, Eq. [§]) can 
be expressed as an integral over the probability density of the independent other particles. 
With the new definition of the interaction parameter, this forces reads 

N 

F * = X^ f j< MkU^iO- (25) 
1=1 



AV Z r r2n 

N z =^fN = n 
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Expressed by the probability density of the N interacting particles it becomes 

Fi = n drj^Tji /idrl^v ji))^. p(n + r jh t) . (26) 
Jv 

Finally, we derive an expression of the projection of the relative velocity. After a few 
calculations, one obtains the relative velocity as a function of the velocity angles tp i} ipj and 
the position angle otf 

/ Pi + Pi \ f Pi ~ %i 
Vji = Vjifai, tfj, oij) = -2s sin I — a s I sin I 1 , (27) 

which allows a closure of the kinetic description. 

By virtue of the assumptions made, the effective force separates in three parts: one 
coming from the repulsion at short distances, one from "moving away" and the third, from 
"approaching" individuals. In the following discussion, we will first restrict ourselves to the 
second and the third interaction type, which corresponds to a limit of a vanishing short-range 
repulsion length l r = 0. In general, the short-ranged repulsion increases the pressure and 
the effective temperature in the system but does not induce instabilities and inhomogeneous 
steady states as observed in the simulations (see Sect HVj) . However, it determines the size 
of the observed structures. 

In the following, we investigate the force turning the velocity angle and project on the 
unit vector e Vi (t). The integral ( 1261) transforms into 

-<»m Jrji<l a 

~ TT / dr ji (ViiO(-Vji)*3i)<Pi P( T i + Y f" l ) } ( 28 ) 

where the spatial distance between the j'-th and the i-th particle appears explicitly in the 
argument of the PDF. The product between the two unit vectors in ( 128]) can be expressed in 
terms of the different angles as e 9 -iji = sin(ay — (pi). Further on, we will also not distinguish 
from the particle number approaching and moving away by setting N m = N a . 

From Eq. (I27j) we find the two different regions of integration where the relative velocity 
has a different sign for tpj in the interval ip j < tpj < 2tt + <pf. 

v H > for fi±£i < a < 2i±£± + n "moving away" , 

(29) 

< for + n <a i^ + 27r "approaching" , 

Hence, a particle, located in the half-sphere in clockwise direction from the mean angle 
(<fj + <fi)/2 approaches the focal particle, whereas, a particle located in the other half-sphere 
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FIG. 3. Visualization of the support in the spatial integration for "approach" and "moving away" 
in the kinetic description for a binary interaction of the i-th individual (particle), moving with 
velocity Vj (heading angle ipi), with a neighbor j within its sensory range l s moving with velocity 
Vj (heading angle (fj). The dotted line, determined by the mean angle {tpi + (fj)/2, represents the 
border between the two distinct spatial regions (half-spheres) corresponding to "approach" and 
"moving away" of individual j: If the relative position vector r^j of individual j with angle ap- 
points above the dotted line, into the red half-sphere, then the two particles are coming closer 
("approach"). If rji points below the dotted line (blue half-sphere) then the two particles move 
away from each other (as in this example). The dashed vectors correspond to the center-of mass 
of the correspond ing half-spheres. 

counterclockwise from the mean angle are moving away. Thus, the support of the interaction 
integrals for "approach" and "moving away" corresponds to these two different half-spheres 
(see Fig [3]). Please note that for (pj = {pi the social force vanishes as Vji = f[2"T|) . 

As an approximation, we will replace this integration by averaged values fixing the prob- 
ability distribution in the centers of the half-spheres, i.e. at a distance from l s /2 and 
with the angle ctj A = (<Pj + V 9 «)/2 — tc/2 for the approaching particles and correspondingly 
a j,m = {<£j + <Pi)/2 + 7r /2 for particles moving away. The corresponding distance vectors are 
parallel but point in a different direction. They read 



r 



Til 



— — e <Pj+<Pi+n , 

Z 2 



S 



a 




(30) 



2 
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These vectors correspond to the center of mass of approaching (moving away) individuals 
under the assumption of a homogeneous spatial distribution of neighbors with direction tpj. 

We separate the spatial PDF from the distribution of the velocity-angles using the con- 
ditional angle PDF, i.e. P(r,tp,t) = p(r,t)P(tp,t\r). Assuming that the density depends 
only weakly on the position within the sensing region, one can take it constant within the 
range of spatial integration. Consequently, the spatial integration cancels together with the 
denominator and we get for the projected force 

ftpi+2-K 



F 4 e w = + s Hm / dipjSmfa - ifi) P(tpj,t\Ti + r m ) (31) 

-s jj, a / dtpj sm(ipj - ifi) P(ipj, t\r,i + r a ) 

The resulting expression will be taken in a dipole approximation. Subsequently, we develop 
the difference in small l s up to the first derivative which yields two contributions F^e^ « 
Fje^ + Fje^. In first order, we get 

F 4 eW = (n m - (i a ) / dtfj sin(ifj - tpi) P(ipj, *|r 4 ) (32) 
= (fl m - fl a ) {co8(<pi)s(Ti,t)-sm((pi)c(Ti,t)) 

The second approximation includes the derivatives in direction of the centers of the half- 
spheres. Inserting the corresponding directional cosines we obtain in second approximation 



F l eg ) = |(/x m + fi a ) 



d d 
cos(<^) ( r m — sin(v9 i )) j - sin(^) ( f m — cosfaj)}. 



(33) 



and the unit vector f m of ( )30|) still depends of the velocity angles which have to be taken 
into account if averaging = J d<pj...P(ipj,t\ri). 

Both orders depend qualitatively different of the social force strengths. Whereas the first 
order, corresponding to a monopole contribution, depends on the difference of the interaction 
strengths, the second (dipole contribution) contains the sum of both coefficients. Assuming 
that both coefficients differ only in their sign, the first order describes effectively a perfectly 

(2\ 

symmetric "escape and pursuit" (or "Head on Head") situation. In this case F^e^ , which 
depends on the local gradients, vanishes. Thus in this special cases, the effective turning 
rate does not depend, up to second order, on spatial (density) inhomogeneities. On the other 
hand for social interactions with equal sign, such as "pure attraction" and "pure repulsion" , 
the F;e^ vanishes and only the second order contributes. In this case, there is no monopole 

13 



force driving the local collective dynamics. However, depending on the sign, the social 
force can lead to an amplification of density inhomogeneities ("pure attraction"), or will 
stabilize the spatially homogeneous state ("pure repulsion"). We should mention that the 
short ranged repulsion is described by a similar expression as for the case of pure repulsion 
but without the factor s . 

In the general case of arbitrary fi a and /i m , the resulting dynamics will be a combination of 
both contributions. However, from the analysis of the dependence of the effective turning as a 
function of the interaction strengths, it can be easily seen that that the main diagonals fi m = 
fx a and n m = —fi a separate qualitatively different regions. For example, in the escape and 
pursuit quadrant, the second order contribution changes its sign at the diagonal with \n a \ = 
\fj, m \, which indicates that the social interactions tend to amplify density inhomogeneities 
in the pursuit dominated region \fj, m \ > \fi a \, whereas the opposite is the case in the escape 
dominated region |/i m | < \fi a \. 

Further on, we will focus for simplicity on the monopole term. The averaged force terms 
in Eqs. (|T3|) become 



(Fie^sincpi) = s (fi m - Li a ) p(r h t) {(T Xuy .(r h t)) s(r h t) - T yuy .(r h t) c(r h t)} , (34) 



(F*ejy cos = s (fj, m - Li a )p(ri,t) {{T XijXi {r h t)) s(r h t) - T Xityi (T U t) c(r h t)} , (35) 

respectively. With these expressions, the term in the evolution equation of the order param- 
eter (Eq. |2T|) containing the forces becomes proportional to (/i. m — fi a ) and we obtain 



(Pm - fJ-a) {c 2 (ri,t)T XitXi (ri,t) + s 2 (ri,t)T yuyi (ri,t) - 2 cfo, t) s(r i; t) T x . >y .(ri, t)} . 



One sees that this term is determined by the temperature of the particle gas. The r.h.s. of 
the given expression can become positive for positive /i m — ii a . It describes the creation of 
order if its value is larger then the effective relaxation process which arises from the noise 
in the velocity angles. 

For an estimation of the critical values of the onset of order one might take the approxi- 
mation in Eq. ( TX5|) which yields that the r.h.s. approaches (/i m — Li a ) (c 2 (rj, t) + s 2 (r i; t))/2. 



and 




(36) 
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Hence, a necessary condition for a growing order parameter, thus onset of collective motion, 
we obtain 

/im - fJ>a > 2^ . (37) 

In general, for vanishing fluctuations, collective motion can emerge only above the main 
diagonal with p m > p a . With increasing noise strength D v the regime of collective motion, 
is predicted to recede into the escape-pursuit regime. 

The above result, as well as the formulation of a kinetic description in general, not only 
provide qualitative insights into the impact of the interaction strengths on the large scale 
system dynamics, but provide also a starting point for a more quantitative analysis of the 
stability of the inhomogeneous solutions. However, this requires further assumptions on the 
properties of the involved probability densities which ensure a closure of the descriptions, 
e.g. assumptions for the temperature, which are beyond the scope of this work. 



IV. SIMULATION RESULTS 



In order to characterize the behavior of the model, we have performed systematic numer- 
ical simulations for varying interactions strengths (—5 < p a<m < 5), different densities p and 
noise strengths D p . 

In the following we will discuss our results in term of the dimensionless density p s = 
Nl^/L 2 , rescaled by the sensory range, which is proportional to the average number of 
individuals per interaction zone for a homogeneous (random) spatial distribution. 

We focus in particular on the question what combinations of p a , p m lead to large scale 
collective motion. The degree of collective motion after the system reaches a steady state is 
measured using the time averaged center-of-mass speed normalized by the preferred speed of 
individuals Sq, which is the well-known order parameter used in the analysis of Vicsek-type 
models: 

(S) t = -{\{vM)t > (38) 
So 

where (-) t denotes the temporal average and the ensemble average. In addition, we 
measure the spatial inhomogeneity (clustering) by the time averaged scaled neighbor number 

{N)t = / (N r (t)+N m (t)+N m (t)) N \ ^ 

\ -^*max / t 
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with (N r (t) + N m (t) + N m (t))N being the average number of neighbors within the metric 
distance given by the sensory range l s of an individual at a given time t. The density- 
dependent scaling number iV max defines the maximal expectation values for the measured 
neighbor numbers corresponding to the closest packing of individuals by assuming an im- 
penetrable repulsion zone with a diameter l r : 

4Z 2 

N max = r }2d -±-l. (40) 

Here r]2d = 7r/(2\/3) ~ 0.907 is the packing fraction for the closest packing of discs in two 
spatial dimensions. The term —1 in the definition of iV max takes into account that the focal 
particle is not being counted as its own neighbor. Please note, that as we are considering a 
soft core interaction, (N) t can in principle be larger than one, in particular for high densities 
and strong attraction. 

Throughout this work we set so — 1, l r — 1 and l s — 5. Furthermore, we use fx r = 20, 
which ensures that for binary interactions the short-ranged repulsion is always larger than 
the sum of the other possibly attracting forces. The particle number is set constant to 
N = 2000 and the density is varied by changing the system size L. The parameter-space 
diagrams in Figs. |4] and |5] were obtained from interpolating the results for (S) t and {N) t 
for 441 (21 x 21) individual, evenly spaced, grid points in the interaction parameter space 
with —5.0 < fi m , fx a < +5.0. Each such point corresponds to an average over the results of 6 
independent simulation runs, whereby for each run a temporal average was taken after the 
system reached a steady-state. 

At sufficiently high densities and sufficiently low (angular) noise, we can observe the 
onset of collective motion for a wide range of interaction parameters (see Figs. H]and[5]). At 
low noise and high density, the region of collective motion coincides approximately with the 
escape and pursuit quadrant of the interaction parameter space with \x m > and fx a < 0. It 
contains the special cases of only attraction to moving away (repulsion from approaching) 
individuals with Li m > and /x a — (ji a < and /x m > 0), and extends also into the 
pure repulsion region (jx a , fi m < and fx a LX m ) and to a much lesser extent into the pure 
attraction region (/i a ,/i m > and fi m fi a ). This is agreement with the predictions of the 
kinetic theory for the monopole approximation the effective social interaction. However, the 
region of collective motion is smaller in simulations than in the simple theory. This may 
be due to the impact of density inhomogeneities, finite short-range repulsion, and/or higher 

1(3 



p s = 0.56, D v = 0.1 

center of mass speed (S)t neighbor number {N) t 




Ha Ha 



FIG. 4. Steady-state center-of-mass speed (S)t (left) and the neighbor number (N) t (right) versus 
fj, m and n a for different densities p s = 0.56 (top), 1.25, 5.00 and = 0.1. The vertical and horizon- 
tal dashed lines indicate the zero axes. The diagonal dashed-line in the escape and pursuit quadrant 
indicates the border between escape-dominated (below the diagonal) and pursuit-dominated be- 
havior (above the diagonal). 
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p a = 1.25, D v = 0.02 

center of mass speed (S)t neighbor number {N) t 




center of mass speed {S)t neighbor number {N)t 




Ha fJ> a 



FIG. 5. Steady-state center-of-mass speed {S)t (left) and the neighbor number (N)t (right) versus 
fj, m and fi a for low angular noise D v = 0.02 (top) and high angular noise = 0.5 (bottom). The 
dashed lines as in Fig. 01 

order effects. 

Within the escape and pursuit regime, where the repulsion to approaching individuals 
(escape) dominates over the attraction to moving away individuals (pursuit), the neighbor 
number is low. This corresponds to low degree of clustering and a rather homogeneous 
spatial distribution of particles throughout the system. The neighbor number (n) increases 
strongly in the pursuit dominated regime (\(J, m \ > \(J, a \ with /i m > 0, fi a < 0) indicating strong 
density inhomogeneities corresponding to dense collectively moving bands and clusters. This 
resembles the behavior observed in the original Brownian particle escape &: pursuit model. 

The neighbor number is also high in the pure attraction regime (/U a ,/i m > 0) without 
collective motion, where clusters with vanishing center of mass velocity can be observed. 
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Interestingly, at moderate densities (p s = 0.56, 1.25 in Figs. I4f5l) . the maximum of the 
neighbor number is located in the pursuit-dominated regime with collective motion, and 
not, as one might expect, in the regime of (strong) overall attraction (fi m , fi a ^> 0). In the 
ordered state, particles move approximately in the same direction and the relative speed \vji\ 
will be close to zero. For low repulsion from approaching individuals the escape response is 
negligible, whereas the attraction to moving away individual suffices to maintain cohesion, in 
particular at low noise strengths. Effectively, the density of such collectively moving cluster 
is limited by the short-range repulsion (see Fig. [6]B). In the pure attraction regime, particles 
on the boundaries of a cluster will be attracted towards the local center of mass. However, 
due to the self-propelled motion with inertia and scattering with other individuals within 
the disordered cluster they will eventually move outwards again. As a result we observe 
disordered aggregates, which resemble mosquito swarms (see Fig. [2]B), and are more dilute 
in comparison to the coherently moving clusters in the pursuit-dominated case. 

As might be expected, increased stochasticity in the motion of individuals, inhibits the 
onset of collective motion. The region of parameter space with (S)t significantly larger than 
reduces strongly with increasing D v by receding towards the regime of strong escape and 
pursuit response (see Fig [5]) in agreement with the kinetic theory. 

For the Head-on-Head regime as well as for pure repulsion (with \i m ^> 0) a quasi homo- 
geneous distribution of particles can be observed with no collective motion (see Figs. [2p,D 
and EH)]). 



V. DISCUSSION 

In this work, we have analyzed a model for collective dynamics based on selective attrac- 
tion and repulsion interactions, which was recently used to model the evolution of phenotypic 



phase change in locusts [32j. The modelling of individual dynamics in terms stochastic dif- 
ferential equation (Langevin equations), allows a straight forward derivation of a kinetic 
description, which may be used for further theoretical analysis based on mean-field consid- 



erations and moment expansion of the corresponding probability density function [35|, |36 |. 

The model is able to account for three types of social responses relevant in the biological 
context: escape and pursuit, pure avoidance and pure attraction behavior. We have shown 
that large-scale collective motion of self-propelled particles emerge without any explicit 
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FIG. 6. Examples of spatial configurations for the pure escape case (A) with \i a = —3 and /x m = 
and the pure pursuit case (B) with/i a = and /i m = +3. 



velocity-alignment mechanism in the (generalized) escape and pursuit regime, without the 



Hi- 



spatial anisotropy in the social interaction in the original escape and pursuit model 
Here, we should note that in the original escape-pursuit model, with Brownian dynamics of 
individual agents, the spatial anisotropy is essential for the emergence of directed collective 
motion. 

In general, the spatial distribution of individuals during collective motion depends 
strongly on the relative strength of the different social forces. If escape dominates we 
observe a homogeneous spatial distribution (Fig. EK), whereas in the pursuit-dominated 
case compact, coherently moving structures, as for example snake-like clusters, can be ob- 
served (Fig. [6j3) . In between, for comparable escape and pursuit strengths the band like 
structures perpendicular to the average direction of motion emerge, which appear also in 
systems with velocity-alignment (Fig. |2]A.) 



The region of collective motion decreases with increasing noise as well as with decreasing 
density. However, at low densities, the region of collective motion shows a clear shift towards 
the pursuit-dominated regime (see Fig H] top), where in a finite system a order parameter is 
maintained by relatively few moving clusters containing most of the individuals. 

These results are not only in agreement with the basic, qualitative predictions, drawn 
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from the kinetic theor y in Sect. 1111} but resemble also the qualitative behavior of the original 



escape-pursuit model 



11] . Interestingly, at intermediate densities and within the interaction 



range studied, the maximum of the neighbor number - indicating strongest clustering in the 
system - appears in the pursuit-dominated regime and not for pure attraction. This can be 
understood from the fact that local order decreases the effective "temperature" associated 
with absolute deviations of the velocities of single particles from the average velocity of 



their neighbors 
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351 ] . As a result we observe a decrease in the active pressure due to the 



stochastic self-propelled nature of individual motion, which counteracts the concentration 
of individuals due to attractive forces. 

In conclusion, the modelling of collective motion in biology via selective attraction- 
repulsion interactions appears very promising. The model accounts for various individual 
behaviors and displays different spatial patterns of collective motion. The response based 
only on the distinction between approaching and moving away individuals can be directly 



linked to the response to 



role in various species 



ooming visual stimuli, which has been shown to play an important 
Finally, in this context, we should mention a recent work by 
They introduce a model for collective motion based on selec- 



37M39) 
rs 4(J- 



Lemasson and coworkers 
tive interaction of individuals based on explicit, simplified description of visual information 
available to each individual. 
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